Computational Approach to Gravitational Waves 
Forms in Stellar Systems as Complex Structures 
through Keplerian Parameters 
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Abstract — In this paper we investigate the gravitational waves 
emission by stellar dynamical structures as complex systems 
in the quadrupole approximation considering bounded and 
unbounded orbits. Precisely, after deriving analytical expressions 
for the gravitational wave luminosity, the total energy output and 
gravitational radiation amplitude, we present a computational 
approach to evaluate the gravitational wave-forms from ellipti- 
cal, circular, parabolic and hyperbolic orbits as a function of 
Keplerian parameters. 

Index Terms — gravitational waves; theory of orbits; numerical 
gravity; 

I. Introduction 

In recent years, detailed information has been achieved for 
kinematics and dynamics of stars moving in the gravitational 
field of such a central object. The statistical properties of 
spatial and kinematical distributions are of particular interest. 
Using them, it is possible to give a quite accurate estimate of 
the mass and the size of the central object. More precisely, in 
ill], it is described a campaign of observations where velocity 
measurements are extremely accurate. Then from this bulk 
of data, considering a field of resolved stars whose proper 
motions are accurately known, one can classify orbital motions 
and deduce, in principle, the rate of production of GWs 
according to the different types of orbits. This motivates this 
paper in which, by a classification of orbits in accordance 
with the conditions of motion, we want to calculate the 
GW luminosity for the different types of stellar encounters. 
Following the method outlined in fT\, (3], f4|, we investigate 
the GW emission by stellar dynamical structures as complex 
systems in the quadrupole approximation considering bounded 
(circular or elliptical) and unbounded (parabolic or hyperbolic) 
orbits. 

The main parameter is the approaching energy of the stars 
in the system (see also [Sl and references therein). However, 
in a recent work 0, it is shown that the parametric form of 
the main evolution equation of Lemaltre-Tolman-Bondi (LTB) 
cosmological model can be obtained considering the limiting 



rectilinear solutions of Keplers problem. More precisely, the 
analogy between the relativistic evolution of LTB shells and 
the classical discussion of a falling body in a Newtonian center 
of attraction, allows us to conclude that the relativistic orbits 
are clearly related to suitable astronomic coordinates. Here, 
these coordinates can be adopted in order to have a detailed 
investigation of the gravitational waves emission by stellar 
encounters. 

Precisely, briefly recalling the main features of stellar en- 
counters and orbit classification, we present the computational 
aspect to obtain the gravitational wave-forms for elhptical, 
parabolic and hyperbolic encounters. 

II. Orbit in Stellar Encounters 

Let us take into account the Newtonian theory of orbits since 
stellar systems, also if at high densities and constituted by 
compact objects, can be usually assumed in Newtonian regime. 
We give here a self-contained summary of the well-known 
orbital types in order to achieve below a clear classification 
of the possible GW emissions. We refer to the text books [71, 
fSl, f9l for a detailed discussion. 

A mass nii is moving in the gravitational potential $ gener- 
ated by a second mass TO2- The vector radius and the polar 
angle depend on time as a consequence of the star motion, i.e. 
r = r(t) and — 4'{t). As it is known, the shape of the orbit 
depends by the total energy H of the system. Precisely 

when H < the orbit is an ellipse and the equation of the 
trajectory is 

^ = T— ^ J ' 

1 + e cos (p 

where I is the so-called semi-latus rectum or the parameter of 
the ellipse and e is the eccentricity of the ellipse (for more 
details see f?]). 

In a generic elliptic Keplerian motion, called a the major 
semiaxis and E the eccentric anomaly, the orbit can be written 
also as (see g), ||9l) 



a{l — ecosE) 



(2) 



hence, there is the following relation between the eccentric 
anomaly and the angle 0: 

cos E — e 



coscf) = -. (3) 

1 — e cos E 

While _ff > is the condition to obtain unbounded orbits. 
The trajectory is 

^ = T^-^ 2 ' 

1 + e cos cf) 

where e > 1. The equal sign corresponds to iJ = 0. Therefore, 
in order to ensure positivity of r, the polar angle has to be 
restricted to the range given by 

1 + ecos0 > 0. 

This means cos0 > 1, i.e. G (— tt, tt) and the trajectory is 
not closed any more. For (f) ±7r, we have r —* oo. The 
curve (|4]i, with e = 1, is a parabola. For e > 1, the allowed 
interval of polar angles is smaller than (p G (— tt, tt), and the 
trajectory is a hyperbola. Such trajectories correspond to non- 
returning objects (for more details see [^). 
Lets consider a generic Hyperbolic Keplerian Orbit, with 
eccentricity e, real semi-axis a and F as variable, analogous 
to the elliptic eccentric anomaly E. The orbit is defined by 

r Q;(ecoshF — 1); (5) 

hence, there is the following relation between F and the angle 

ct>: 

^ ^ — Q;(ecoshF — 1) 
ea{e coshF — 1) 

Lets consider the Keplerian Parabolic Orbit. It is defined by 
the relation 

p2 

where P is a parameter In this case 

21 - P2 



(7) 



p2 



(8) 



In the next section we present the computational 
approach to obtain gravitational wave-forms in the different 
cases as function of the eccentric anomaly, the angle F and 
the parameter P. 

Remark 1. Let the eccentricity tend to unity in equations 
(|7); in the limit the corresponding orbits are called rectilinear 
ellipse, hyperbola and parabola respectively. 
The evolution of the r-shell in LTB models is analogous to 
a Keplerian motion on a rectilinear ellipse, hyperbola or 
parabola (for more details see 



III. Computational Approach to the Gravitational 
Wave Forms 

At this point, considering the orbit equations, we want 
to classify the gravitational radiation for the different stellar 
encounters. We send the Reader to the References ifTOl . ifTTI . 
ifTZI . lfT3l for a detailed exposition. 

The Einstein field equations give a description of how the 
curvature of space-time is related to the energy-momentum 



distribution. In the weak field approximation, moving massive 
objects produce gravitational waves which propagate in the 
vacuum with the speed of light. One can search for wave 
solutions of field equation from a system of masses undergoing 
arbitrary motions, and then obtain the power radiated [4). 
The result, assuming the source dimensions very small with 
respect to the wavelengths (quadrupole approximation |[8]), is 
that the power ^ radiated in a solid angle Vl with tensor of 
polarization is 

dH _ G (d^Q.j 



Bttc^ \ dt^ 



where is the quadrupole mass tensor 



G being the Newton constant, ra the modulus of the vector 
radius of the a-th particle and the sum running over all masses 
TTia in the system. With this formalism, it is possible to 
estimate the amount of energy emitted in the form of GWs 
from a system of massive objects interacting among them f2|, 
|3 1. In this case, the components of the quadrupole mass tensor 
in the equatorial plane [0 — ^) are: 



ixx 



pr'^{3 cos^0 - 1) , 



]yy = ^r^(3sin^0 - 1) 

„2 



; zz 



-fir 



(9) 



ixz ^ zx 



Qyz — Qz 



, 



, 



Qxy = Qyx = ipr'^ COS (j) s'm (/> , 

where n ~ '"^"^ is the reduced mass of the sys- 
tem and the masses mi and rn2 have polar coordinates 
{r,; cos 9 cos (f), Ti COS Osincj), sin 9} with i ~ 1,2. The origin 
of the motions is taken at the center of mass. 

Direct signatures of gravitational radiation are its amplitude 
and its wave-form. In other words, the identification of a GW 
signal is strictly related to the accurate selection of the shape of 
wave-forms by interferometers or any possible detection tool. 
Such an achievement could give information on the nature of 
the GW source, on the propagating medium, and , in principle, 
on the gravitational theory producing such a radiation |14|. It 
is well known that the amplitude h of GWs can be evaluated 
by 

on 

hjk{t,R)^j^Qjk, (10) 

R being the distance between the source and the observer and 
{j,fc} = l,2. 

Moreover, in accordance to 



{hi. 



^22 



2hl,) 



1/2 



(11) 



we can compute the expected strain amplitude for different 
orbits. 

Note that we are going to study the evolution of compact 
binary systems that are formed through the capture of a 
moving mass nii by a gravitational field, whose source is a 
Black Hole(BK) of mass m2 and being mi <^ m2- So we 
will have that the reduced mass of the newborn binary system 
II ~ mi and the mass ratio — is ^ 1. This constraint is 

^ ^ "1-2 

satisfied by several real systems. 

Let us now derive the GW amplitude in relation to the 
orbital shape of the binary systems. 

A. Quadrupole and GW Amplitude from Elliptical Orbits 

In this section the components of the quadrupole mass 
tensor in the equatorial plane for the elliptical orbits are 
expressed. Precisely, by Eq.s (|2|i, (O and (|9|l, we obtain: 

Q^^ = -a'^m.i [-Se^ + 4e cos E + {e^ - 3) cos^ -B + l] , 




— £=0.7 

- -£=0.5 
- e=0.2 



Waveform tor an elliptic orbit 



Fig. 1 . The gravitational wave-forms from elliptical orbits shown as function 
of E{t). We have fixed the masses of the order IAMq. m2 is considered at 
rest while mi is moving. The distance of the GW source is assumed to be 
R = 8 kpc and the eccentricity is e = 0.2, 0.5, 0.7. 



Qyy = o?mx [-3£2 + 2e cos E + {le^ - 3) cos^ + 2] 
Qzz = — a^mi(ecos£' — 1)^ , 

Qxy — Qyx — 



= 3a'^mi{e - cos £;) (e cos - ^) J ~ co1^e-i)^ 



(12) 

Moreover, considering a binary system and the single compo- 
nents of Eq. ( [Tol l, it is straightforward to show that 

_ c(Givm2nii[(e-l)(e+3)csc^(f )+2(e^-3)cos£;] 
"■11 ~ cosB-1 ' 



l22 



aGivm2mi[(2e^+e-3) csc^ (f )-|-(4e^ -6) cos E] 



iE-1 



(13) 



1 , ory / (e^~l)sin^E 

hi2 - h2i = 3GNm2miaJ- '^Jj^_^., 



Wavetorm for hypertiolic orbit 




Fig. 2. The gravitational wave-forms for hyperbolic encounters as function 
of the F(t). As above, we have fixed the masses of the order IAMq. m2 is 
considered at rest while mi is moving. The distance of the source is assumed 
at ii = 8 kpc. The eccentricity is assumed with the values £ = 1.2, 1.5, 1.7 



(ecos£;-l)(e-4cos£;+cos2£;+2) cot^(.|) csc^ E 

cosB-1 ' 

and then the expected strain amplitude h, which strictly 
depends on the initial conditions of the stellar encounter, can 
be evaluated in accordance to ( fTTT l. 

B. Quadrupole and GW Amplitude from Hyperbolic Orbits 

When hyperbolic orbits are taken into account, Eq.s (|5]l, (|6]l 
and ^ become: 



a^TOi(ecoshi^ — 1) 



3{l+a-aecoshF) 
(ae^ cosh F—ae)'^ 



1 



Qyy — a^mi{ecoshF — 1)^ 



3 1 



(i+g-aecoshF)" 
(ae^ cosh i^ — ae)"^ 



In this case the single components of Eq. ( fTOl ) for a 
hyperbolic orbit, are: 



hu 



2Gjv^Ti2^^ 
7^ 



3;^ gr 

a^(ecoshF-l)^ Q(e cosh(F(t))-l)^ 

i^+a(-2s^ + 8s^-3 



+ 



2Gjv?7i2^i 
7^ 



3;^ _i 91^ 

a-'{e cosh F a{e coshF-1)^ 

ig^+.(4e^-10s^ + 3) 



Qz 



-a^^(ecoshF — 1)^ 



Qxy — Qyx 



3afj.{<i cosh F-l){l+a-ae cosh F)^jl~ ^' + ° 



(14) 



hi2 = h2i 



f (i + Q-oe cosh F)2 . ,\ 
ae{s cosh F — l)'^ 



3/2 



P + a^ (e^ - 1) +a^e (e^ - l) coshF(£ cosh(F - 2)) , 

(15) 

which, as before, strictly depends on the initial conditions of 
the stellar encounter. 




Fig. 3. The gravitational wave-forms for a parabolic encounter as a function 
of P{t). As above, we have fixed the masses of the order 1.4A/q and m2 
is considered at rest, while mi is moving. The distance of the GW source is 
assumed at R = 8 kpc. The eccentricity is e = 1. 



C. Quadrupole and GW Amplitude from Parabolic Orbits 

Finally, considering parabolic orbits, the explicit expression 
of the components of the quadrupole mass tensor in the 
equatorial plane through the parameter P are deduced by (|7|l, 
® and ©: 



Qxx — 



p3— 



1 



p4 



_ lp4 f 3[P'-i2l-l f] A 

'^yy ~ ^ \ pW r ™i 



4 



(16) 



Qxz — Qz 








comparable to the Chandrasekhar limit I.4A/0). This 
choice is motivated by the fact that ground-based experiments 
like VIRGO or LIGO expect to detect typical GW emissions 
from the dynamics of these objects or from binary systems 
composed by them (see e.g. [12]). We would like to underline 
that the presented approach provides the necessary support for 
the application of computational schemes to obtain qualitative 
description of the Gravitational Waves. 
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P'i-(2/-l)2 

pi — —mi 



Qxy — Qyx — 4 (^^ 1)-^^ 

Moreover, the single components of Eq. (fTOl l for a parabolic 
orbit, are: 

L _ iGN'm2mi 
fill — p5 ) 



i22 



8Gjvm2mi 



p2 



h21 



6Gm2mi[(l-2/)P''-(2Z-l)^] 



pa 



(17) 



and then the expected strain amplitude h, which strictly 
depends on the initial conditions of the stellar encounter, can 
be evaluated in accordance to (fTTT l. 

IV. Conclusion 

We have analyzed the gravitational wave emission coming 
from stellar encounters in Newtonian regime and in quadrupole 
approximation. In particular, we have taken into account the 
expected strain amplitude of gravitational radiation produced 
in tight impacts where two compact objects with masses 



